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SUMMARY 


The low temperature provided by 3 He refrigerators (0.3-3 K) have useful space 
applications. However, the low temperatures and the low surface tension of 3 He 
require special design considerations. These considerations include the need for 
small pores to contain the liquid in a matrix, the effects of bubble rucleation 
and growth, and the effects of the thermal conductivity within the matrix. These 
design considerations are discussed here along with an analysis of a possible 
confinement system. 


SYMBOLS 


a acceleration 

B Bond number 

o 

B x B /rh 

o 

E election number 

J 

Ex E^r/Qh 

F p Va: force due to acceleration 

a a 

F force due to a pressure difference 

P 

F force due to surface tension 

s 

F force due to surface tension gradients 

g gravitational acceleration 

h height of liquid 

h Planck constant 

k Boltzmann constant 

K effective conductance 

e 

conductance of liquid 
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K conductance of matrix 

m 

L latent heat 

n number of capillaries 

N number of liquid molecules 

P pressure 

P F /-nr 2 ; pressure due to surface tension 

s s 

AP pressure difference 

Q heat flux 

Q total heat flux 

P 

r radius 

R gas constant 

t hold time 

T temperature 

condensation temperature 
operating temperature 
superheat 

T saturation temperature 

sat 

AT temperature difference 

V vo lume 

x distance m the direction of a 

a rh 

3 Qh/r 

e mass fraction lost during pump-down 

k effective conductivity 

e 

k ^ liquid conductivity 

k matrix conductivity 

m 

n volume void fraction 

liquid density 
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P v vapor density 

a surface tension 


INTRODUCTION 


The development of a space-compatible 3 He refrigerator would provide a signifi- 
cant improvement m several areas of research such as IR astronomy and the study of 
critical phenomena in quantum fluids (refs. 1 and 2). 3 He can produce refrigeration 
by evaporative cooling from 3.2 K (its critical point) to about 0.3 K. The lower 
limit is set by the vapor pressure and the system pumping speed. 

While He refrigeration has been used for a number of years in laboratories 
(refs. 3 and 4), on balloons (ref. 5), m aircraft (ref. 6), and in spin-stabilized 
sounding rockets (ref. 7), it has yet to be used in the low-gravity environment of 
space. Unlike ‘‘He, which was used in IRAS (ref. 8), 3 He is not a superfluid at 
these temperatures and cannot be contained by the fountain pressure generated by a 
temperature gradient. 

A proposed 3 He refrigerator is shown in figure 1. While this is not the only 
possible configuration, this design will be used here as a reference m discussing 
the major features that are also common to other designs. This design consists of 
an adsorption pump, two heat switches, and a pot. The pot, the pump, and the inter- 
connecting tubing form a sealed system that need be filled with 3 He only once. The 
heat switches connect the refrigerator to heat sinks whose temperatures are indi- 
cated in the figure. The refrigeration cycle starts with the system below 10 K. 

The pot switch is closed and the pump switch is opened. The adsorption pump is then 
heated to drive out the 3 He gas which condenses in the pot. When the condensation 
is complete, the heater is turned off and the switches are reversed. The pump 
cools, then readsorbs the gas, causing evaporative cooling m the pot. When the 
pot runs dry, the cycle can be repeated. 

Variations of this type of refrigerator have been used in the past whenever 
gravity or centrifugal force ensures that the liquid remains in the pot. The pur- 
pose of this paper is to discuss various aspects of containing the 3 He m space. 

If the pot is filled with a porous material the liquid 3 He can be contained by 
capillary attraction. Although surface tension methods have been used for control- 
ling other fluids m space, the low temperatures and low surface tension of 3 He 
require special considerations. The use of capillary confinement in a 3 He refriger- 
ator was first suggested by Ostermeier (ref. 9). Since then, Ennis (ref. 10) has 
demonstrated condensation, confinement, and useful refrigeration of both 4 He and 
3 He in an inverted geometry. In many ways the inverted operation is a more severe 
test than a zero-gravity operation. 

The main concerns with capillary confinement are whether the vapor can be 
condensed into the porous matrix without leaving large voids, and whether the liquid 
can be evaporated from the sponge without forming bubbles that would expel liquid. 
Expelled liquid would reduce the quantity of liquid available for cooling. 

These concerns are discussed in the following sections, followed by a model 
of a proposed refrigeration system. The porous matrix will be modeled as a cluster 
of smooth cylindrical capillaries. These capillaries will be assumed to be parallel 
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to any acceleration of the refrigerator. Although real accelerations can occur in 
any direction, accelerations parallel to the capillary axis produce the largest 
effect. Deficiencies in this model will be discussed. 

We would like to acknowledge the support of the NASA Office of Aeronautics and 
Space Technology. 


CAPILLARY CONFINEMENT 


The porous matrix serves two functions during the refrigeration cycle. It pro- 
vides a place for the vapor to condense during the condensation phase, and it 
retains the liquid during the evaporation phase. 

In order for the vapor to condense, several conditions must exist. The part of 
the system which is m contact with the vapor must be at a temperature below the 
critical temperature of the gas, and the pressure in the system must be greater than 
the corresponding vapor pressure. If these conditions are met, condensation will 
occur. If several areas are cold enough for condensation to occur, the condensation 
can be expected to occur preferentially at the coldest point. This effect is due to 
the vapor pressure decreasing and the surface tension increasing with decreasing 
temperature. Both the vapor pressure gradient and the surface tension gradient pro- 
vide forces that drive the condensed liquid to the coldest point (cold spot) in the 
system. Thus, at the end of the condensation phase, all of the liquid will be con- 
densed at the coldest point unless there is some other force, such as gravity, to 
drive the fluid elsewhere. 


If the cold spot is a porous matrix, the condensation is enhanced because the 
surface tension is increased by the small size of the pores, and because the 
increased surface area increases the condensation rate. However, there may be a 
difficulty when using a porous matrix — that the condensation might take place only 
on the outside of the matrix, leaving voids in the interior. However, experiments 
by Donnelly (ref. 11) and Ennis (ref. 10) have not demonstrated this effect. They 
showed that the condensed liquid substantially fills the matrix. 


In a spacecraft, the porous matrix must hold the liquid against any lateral 
accelerations. This is required during both the condensation and evaporation phases 
of the refrigeration cycle. These accelerations can, in general, occur in any 
direction and can be of variable magnitude. The effect of the acceleration forces 
is best described in terms of the Bond number, B (refs. 12 and 13). The Bond 
number is the ratio of the acceleration forces to the surface tension forces: 


B 

o 


fa _ V rh 

F 20 

s 


( 1 ) 


where we have assumed that the liquid wets the capillary wall with a contact angle 
of 0°; so the cosine term that normally appears in the expression for the surface 
tension force can be ignored. This assumption is valid for 3 He (ref. 9) and will be 
used throughout this paper. If B >1, the acceleration forces dominate, and if 
B <1, the surface tension forces dominate. Thus, to retain the fluid, the matrix 
must be selected such that the stability condition, B << 1 is met for all 
expected accelerations. Since and a are properties of the refrigerant, and 
a is determined by the environment, the designer is free to choose only r and h. 
This restriction can be emphasized by separating B^ into two factors: 
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( 2 ) 


B = aBi 
o 

where Bx = pa/2 a gives the system properties and a = rh gives the" design param- 
eters. The stability condition can now be written as a << 1/Bi. Since Bi is a 
function of temperature, the stability condition must be evaluated over the whole 
operating temperature range, which is the range from the condensation temperature to 
the evaporation temperature. For the case in which a = 9.8 m/sec (1 g) , the func- 
tion 1/Bi for 3 He is shown m figure 2. For the curve m figure 2, the area 
below the curve is the stable region. For many applications, the accelerations can 
occur in any direction. Thus h must be the maximum dimension of the matrix from 
an open pore. This would lead one to think that the largest possible volume would 
be a sphere of radius h. However, the matrix is not floating free; rather it is 
almost entirely contained. The only opening is in the vent tube (fig. 3) and no 
part of the matrix can be a distance greater than h from the opening of the tube. 
Thus the greatest volume is a sphere of radius h. From the stability condition 
(a << 1/Bi), it is readily seen that the smaller r is, the greater h can be; 
and thus the more liquid that can be held. 


BUBBLE DYNAMICS 


The porous matrix must not only confine the liquid during the condensation 
phase, but also during the evaporation phase. Bubbles that form during evaporation 
tend to expand. At first they will expand radially until they block the capillary. 
If the bubble expands past this point, it will displace and eventually expel 
liquid (ref. 14). Since the expelled liquid will not produce any useful cooling m 
the system, it is desirable to prevent the expulsion. This problem can be broken 
down into several aspects, starting with the initial formation of the bubbles, then 
their growth and movement within the capillaries. This process will be discussed 
m this section along with effects of real nonidealized porous matrices. 


In a porous matrix, with its large irregular surface, bubbles will form by 
inhomogeneous nucleation at nucleation sites on the interface (refs. 15 and 16). 
Bubbles will form when the liquid temperature, T, is raised sufficiently above the 
saturation temperature, T sat (ref. 17). The quantity (T s = T - T sat ) is called 
the superheat. The amount of superheat required for nucleation depends on the sur- 
face properties of the interface. While it is difficult to predict this superheat 
exactly, an estimate can be made. Bald (ref. 16) derived an expression: 


T 

s 



£n < 




16tto 3 

Hn(6N 2 / 3 kT/h) 



(3) 


This function is shown in figure 4 as a function of temperature for 3 He for 1 mole 
of liquid. 


Bubble formation can be prevented if the temperature gradients in the fluid are 
kept small and the superheat is never reached. These limitations require the 
effective thermal conductance K across the fluid to be large. In other words, 
the combined conductances of the matrix and the fluid must be considered between the 
heat source (where nucleation is most likely to occur) and the fluid's free surface 
(where evaporative cooling occurs). If the heat flow is rectilinear, then 
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(4) 


K 

e 


= K + K, 
m i 


where Kju and are the thermal conductances of the matrix and the liquid, 

respectively. AT = Q/K e gives the temperature drop across the system which must be 
less than the superheat. Thus for the system to be stable against the formation of 
bubbles, the following relation must hold: 

Q /K e < T s (5) 


This formula requires that K e be large. Since K.£ is set by the properties of 
3 He, only Kjj, is a free parameter. A high-conductance matrix,- such as copper, will 
give the best stability. The area-to-length ratios in and K.£ are also impor- 
tant. Because the conductivity of copper is considerably higher than 3 He, thick- 
walled capillaries will be more stable than thm-walled ones. 


If we consider the heat flux into a single capillary of radius r and height 

h, then equation (5) can be expressed in terms of an equivalent conductivity < e 

where K = k irr 2 /h: 
e e 


, k iris 
h e 

— << ; 

r 4 - A 


( 6 ) 


(The equivalent conductivity can also be written in terms of the liquid and matrix 
conductivities, and K m , respectively, and the void fraction of the matrix: 

K = < 0 + - n)/n.) 

If a bubble does form, it will be free to move within the liquid. In the 
absence of gravity, which gives rise to buoyancy, the only force available to move a 
bubble comes from surface tension gradients. Such motion is called Marangoni flow 
(ref. 18). The surface tension gradients are the result of temperature gradients 
and the temperature dependence of the surface tension. The surface tension of 3 He 
as a function of temperature is shown in figure 5. Young et al. showed that the 
buoyancy of a bubble will balance the Marangoni effect when 

da dT = 2p £ gr . 

dT dx 3 K } 


Since this balance occurs when the surface tension gradient force, F , is equal to 

F for a = g, F can be deduced: 
a t 


F 


t 


7t r 


2 da dT 
dT dx 


( 8 ) 


Thus the bubble tends to move toward the higher temperature; i.e., toward the heat 
source and away from the free surface of the liquid which is being evaporatively 
cooled . 


Since the bubble most likely forms near the heat source, it will not move far, 
if at all. Furthermore, at low temperatures do/dT approaches zero as does F t . 
Thus at the operating temperature there will be no force to move a bubble. At 
higher temperatures there is another effect that prevents the bubbles from moving. 

At temperatures at which there is a substantial vapor pressure, the bubble acts like 
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a heat pipe (ref. 19). Vapor is evaporated from the hot side of the bubble and is 
condensed on the cold side. This increases the thermal conductivity of the bubble. 
The temperature gradient across the bubble is reduced, thus reducing F . 

We have seen that once formed a bubble is not likely to move; it can only grow. 
Any growth must be accompanied by the displacement and eventual expulsion of liquid. 
This effect can be prevented if the growth of the bubbles can be controlled. To 
illustrate this, we will consider a bubble blocking an otherwise full capillary 
(fig. 6). There are two forces acting on the liquid. One, Fp , is the force due to 
the pressure difference, AP, across the liquid: 


F = rr 2 AP (9) 

P 

The bubble's formation and growth is the result of a heat influx into the bubble. 

This results in a temperature gradient through the liquid. If quasi-static process 
is assumed where the vapor in the bubble and the vapor outside the capillary are m 
thermal equilibrium with their respective liquid interfaces, then equation (9) can be 
written as 


F 

P 


Trr 


dP 

dt 


svp 


AT 


( 10 ) 


The derivative, dP/dT, is evaluated along the saturated vapor pressure curve. This 
can be done in terms of the Clausius-Clapeyron equation (ref. 20): 


dP/dT | sv p = p v L/T 


(ID 


The AT in equation (10) can be expressed in terms of the equivalent thermal con- 
ductance, K : 
e 

AT = Q /K (12) 

where K e = 4- K,, and Q is the heat flux per capillary. As before 

K g = < e rr 2 /h. Using equations (11) and (12), equation (10) can be written as 


F 

P 


Qhp L 
^ v 

K T 


(13) 


The other force acting on the liquid is the surface tension force: 

F = 2irro (14) 

s 

To aid in determining which of the two forces dominates, we will define an ejection 
number, E , as the ratio of F to F : 

1 P s 


E 

1 


Qhp^L 

2rrK oT 
e 


(15) 


If Ej < 1, the surface tension force will dominate and the system will be stable 
against this form of expulsion. The expulsion number may be written as E^ = BEi 
where 
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( 16 ) 


P L 

F = v 

1 2ttok T 
e 


contains the system properties and 


3 = Qh/r 


(17) 


contains the design parameters. The stability condition now becomes 3 < 1/Ei. 

A plot of 1/Ei as a function of temperature is shown m figure 7 for 3 He in a 50% 
copper matrix. This figure shows the regions of stability and instability. 

A real matrix is not likely to be composed of smooth cylindrical capillaries 
but rather of irregular interconnected passages. The effect of the interconnections 
will be to increase the system's stability. The pressure in the bubble will be 
opposed by the surface tension of many pores. The pressure difference across the 
liquid is found by combining equations (9) and (13) : 

Qhp L 

AP = — (18) 

ur 2 K T 
e 


The pressure due to the surface tension at n pores is 


nF „ 

p = s_ = 2no 

s irr 2 r 


(19) 


Thus the Ej = AP/P S is decreased by a factor of 1/n, increasing the stability of 
the system. 


APPLICATION 


In this section we will discuss how the above considerations can be applied to 
the design of a refrigerator. To illustrate, we will consider a hypothetical 
refrigerator for cooling an IR detector. The refrigerator requirements are shown in 
table 1. We will assume that the liquid chamber is a sphere (fig. 3) with a volume 
V = 4Trh 3 /3. Since the capillary device is used both to contain the liquid and to 
conduct heat to the liquid-vapor interface, the entire volume will be filled with 
a high-conductivity matrix such as copper. For our example, we will assume that 
the void fraction is n = 0.5. Since the condensation occurs at T c , a fraction, 
e, of the liquid will be lost in pumping down to T 0 . Furthermore, the volume of 
liquid will change because of the density change between T c and T Q . After pump- 
down, a volume of liquid of Qpt/Lp^ must remain to meet the hold time requirement. 
Thus the volume of the liquid chamber must be 


V 


V 


n(l - e) Lp £ (T c ) 


( 20 ) 


or 
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( 21 ) 


h > J IE 

n -[4irn(l - e) Lp £ (T c ) 

For our example, these are V > 7.8 cm 3 and h > 1.2 cm. 

Bond number considerations lead to another restriction on the chamber specifica- 
tions. This is the requirement that a << 1/Bi, over the entire operating range 
(T c to T q ). We must use the smallest value of 1/Bi over this range. Since 1/Bi 
decreases with increasing temperature (fig. 2), the minimum occurs at T c . Thus we 
have the restriction that 




rh << l/Bi(T ) (22) 

c 

where Bi is evaluated at the peak expected acceleration. In our example the 
requirement becomes rh << 2x10 6 m 2 . 

A third restriction arises from the desire to avoid bubble nucleation. Equa- 
tion (6) gives an expression for h/r m terms of the heat flux per capillary. To 
convert this to an expression in terms of Qp = nQ, we must first estimate the total 
number of capillaries in the system. One such estimate is the ratio of the void 
volume to the volume of a single capillary: 

n = nV/irr z h = 4nh 2 /3r 2 (23) 


Substituting this into equation (6) gives 

3Qp 

^ 4 lynx T 

e s 


(24) 


For our example this gives h > 0.7 mm. Since this is so much less than the pre- 
vious constraint on h (eq. 20), it may be ignored. 

The last restriction comes from consideration of the ejection number. Again we 
must consider the entire temperature range. In this case 1/Ei decreases with 
increasing temperature (fig. 7), so its minimum value occurs at T c . With this con- 
sideration, and again taking into account the number of capillaries in the system, 
the restriction that Qh/r < 1/Ei becomes 

r/h < — (25) 

v/v 

For our example this is r/h < 0.02. 

The restrictions due to hold time. Bond number, and ejection number are shown 
in figure 8. As can be seen in the figure, there is a region of r-h space in 
which a stable refrigerator could be designed. The curve for ejection criterion was 
calculated assuming noninterconnecting capillaries. If interconnecting capillaries 
are used, then this curve would move to the right and down, increasing the stable 
region. The figure shows the case for a single capillary. In practice, n will be 
the number of capillaries exposed to the pump-out tube. 
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CONCLUSION 


It is feasible to design a 3 He refrigerator that will operate in the low-gravity 
environment of space where accelerations can occur in random directions. The key to 
the design is to fill the liquid chamber with a high-conductivity porous material 
such as sintered copper. The pores allow surface tension forces to contain the 
liquid while the conductivity ensures that the vaporization occurs at the surface of 
the matrix rather than internally. The matrix also provides a favorable place for 
condensation to occur and suppresses bubble nucleation and movement. 
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TABLE 1.- REFRIGERATOR REQUIREMENTS 


Parameter 

Symbol 

Value 

Condensation temperature 

T 

c 

2 K 

Operating temperature 

T 

o 

0.3 K 

Hold time 

t 

54 ks 

Refrigeration power Q 

% 

40 pW 

Peak acceleration 

a 

0.1 g 


^Includes parasitic loads. 
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Figure 1.- A zero-gravity 3 He refrigerator showing the principal components. 
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Figure 2.- Plot of 1/Bi as a function of temperature for 3 He for the case in which 
the acceleration is 9.8 m/sec. The region below the curve is the region in which 
surface tension forces dominate. 


14 







Figure 6.- A vapor bubble blocking an otherwise full capillary. The capillary has 
a radius of r. A column of liquid of height h separates the bubble from free 
space. A temperature difference of AT and a pressure difference of AP are 
across the liquid column. There is a heat flux of Q into the bubble. 
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Figure 7.- Plot of 1/Ei as a function of temperature where it has been assumed 
that half of the thermal conductivity path is through copper and half is through 
3 He . 
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Figure 8.- Plot showing the stable and unstable regions of r-h space for various 
types of instabilities. The shading indicates the region of instability. 
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